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Abstract This paper considers contact processes with additional voter model 
dynamics. For such models, results of Lloyd and Sudbury can be applied to 
find a self-duality, as well as dualities and thinning relations with systems 
of random walks with annihilation, branching, coalescence, and deaths. We 
show that similar relations, which are known from the literature for certain 
interacting SDK's, can be derived as local mean field limits of the relations of 
Lloyd and Sudbury. 
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1 Introduction 

Lloyd and Sudbury |SL95[ ISL97L ISudOO| have studied dualities for general spin systems which 
have only two-spin interactions and for which the uniform zero configuration is a trap. Based 
on algebraic considerations coming from quantum theory, they associate dualities with certain 
linear operators. The requirement that these operators are the product of operators acting 
only locally on each site then leads them to consider only duality functions of the form 

M^,y) = ll^''^'^''^'\ (1.1) 

where A is a lattice, x,y G {0,1} are spin configurations, and rj G IR\{1} is a parameter. 
They show that there are lots of dualities between the models they consider. Moreover, they 
show that if two models are dual to the same model (albeit with a different duality parameter), 
then one of these models is a thinning of the other. 

After reviewing some main results of Lloyd and Sudbury in Section [21 in Section |31 we turn 
our attention to contact processes with additional voter model dynamics. Using Lloyd-Sudbury 
theory, we show that such models are self-dual, and moreover dual to systems of random walks 
with annihilation, branching, coalescence, and deaths. We also show that the latter models 
are thinnings of each other, and of the contact-voter models. In Section0]we consider systems 
of interacting SDE's used in population dynamics. More precisely, we consider a version of 
the stepping stone model with selection and mutation, as well as the super random walk with 
an additional quadratic killing. We show that such systems can be derived as 'local mean 
field limits' of contact-voter models, and derive (mostly well-known) dualities, thinnings, and 
Poissonization relations for such models as limits of the Lloyd-Sudbury relations. 

2 Lloyd-Sudbury theory 

2.1 Lloyd-Sudbury dualities 

Two Markov processes X = {Xt)t>o and Y = {Yt)t>o with state spaces Ex and Ey are 
usually -"^ called dual to each other with duality function ip{-,-) if 

E[^P{Xt,Yo)] = E[i^{Xo,Yt)] (t > 0), (2.1) 

whenever X and Y are independent, for arbitrary initial laws P[Xo E • ] and P[Yq S • ] on 
Ex and Ey, respectively. If the functions {ip{-,y) : y S Ey} and {ip{x,-) : y £ Ex} are 
distribution determining, such a duality is informative. 

We will in particular be interested in the case that Ex = Ey = {0, 1} , where A is a finite 
or countably infinite set, and ^p is of the form Hl.l|) . We consider Markov processes in {0, 1} 
which have only two-spin interactions and for which the uniform configuration is a trap. 
More precisely, we consider Markov processes X = {Xt)t>o in {0, 1} with formal generator 

^In [SLQSllSLQ?! . however, the word duahty is used in a much more restricted meaning. 



of the form 

GLsf{x) := Y^ g(i, j)|iax(i)x(j){/(x - 6i - 6j) - f{x)} 

+bx{i)il-x{j)){f{x + 6j)-f{x)} 

+cx{i)x{j){f{x - 6j) - fix)} (2.2) 

+dx{i){l - x{j)){f{x - 6i) - fix)} 



+ex«(l - xij)){fix -6, + 6j) - fix)}], 



where a, b, c, d, e are nonnegative constants and giAxA^Misa nonnegative function such 
that 

Qihj) = qiJ,i) and ^ g(i,j) = 1. (2.3) 

Here 5iij) := 1 if i = j and 6iij) := otherwise, and we adopt the convention that sums or 
suprema over i,j always run over A, unless stated otherwise. The letters a,b,c,d,e denote 
annihilation, branching, coalescence, death, and exclusion, respectively. We may interpret q 
as the jump rates of a continuous time Markov process on A, called the underlying motion of 
the interacting particle system X. The processes in (|2.'2() may be constructed via a graphical 
representation or via their generator; see [Lig85. as a general reference. We call the Markov 
process X with formal generator (|2.2|1 the Lloyd-Sudbury model with underlying motion kernel 
q and parameters a, 6, c, d, e, or shortly, the (g, a, b, c, d, e)-LSM. 

Since there is no spontaneous creation of particles, it is easy to see that a iq, a, b, c, d, e)- 
LSM X started in 

^fin({0, 1}) := {x G {0, 1}^ : ^ x(i) < oo} (2.4) 

i 

stays in this space. We will say that two Lloyd-Sudbury models X and X' are dual to each 
other with duality parameter r] G ]R\{1} if 

E[llr,^t(')^oii)] = E[llr,^oii)Xiii)y (2.5) 

i i 

for all deterministic initial states Xq and X'q. (If X and X' are independent, then integrating 
over the initial laws we see that (|2.5|) holds for nondeterministic initial states as well.) Here 
we assume that Xq,Xq G {0, 1} are arbitrary if A is finite or if |?7| < 1. In case rj = —1 and 
A is infinite, we assume that either Yli-^oii) < cxd or X^j^o(^) < °^- '^^ avoid convergence 
problems, in case |?7| > 1, we assume that |A| < 00. 

We cite the next proposition from |Sudnn| formula (9)] , which is a simplification of |SL95[ 
formula (21)]. While the latter applies only to symmetric models, their formula (20) ap- 
plies to asymmetric models as well. A simplification of that formula in the spirit of |Sudnn| 
formula (9)], plus an outline of its proof, can be found in Appendix El 

Proposition! (Dualities between Lloyd-Sudbury models) The iq,a,b,c,d,e)-LSM 
and iq, a' ,b' , c' , d' , e')-LSM are dual with duality parameter r] £ M\{1}, provided that 

a' = a + 2t/7, b' = b + j, c' = c — (1 -|- t])^, d! = d + ^, e' = e — 7, (2.6) 



where 7 := {a+c— d+rjb) / {1—r]) . In particular, each model with b > andrj := {d—a—c)/b / 1 
is self-dual with parameter rj. 

2.2 Thinnings 

In this section, we recall from |SL97j that if two {q, a, b, c, d, e)-LSM's have the same dual (in 
the sense of Proposition ^), then one is a thinning of the other. Since we will need this further 
on, we define thinnings here in a somewhat greater generality than is needed at this point. 

We may interpret the space N as the space of all particle configurations x on A, where each 
site i G A can be occupied by x{i) = 0, 1, . . . particles. A u-thinning of a particle configuration 
X G N is then obtained by independently throwing away particles, where a particle at i G A 
is kept with probability u{i). More formally, for given x G N and u G [0, 1] , we can choose 
i„ G A such that x = J2n=i ^in (where m is allowed to be 00), and we can choose independent 
{0, l}-valued random variables Xn with P[xn = 1] = u{in). Setting 

m 

Thm^{x):=^XnSi^ (2.7) 

n=l 

then defines a it-thinning of x. We will usually only be interested in the law of Thintj(a;), and 
use this symbol for any random variable with the same law as in 1)2. 7|) . li X,U are random 
variables with values in N and [0, 1] , respectively, then we use the symbol Th.mij{X) for 
any random variable with the law 

P[Thmu{X) G •] := f P[U G du, AT G dx]P[Thin„(x) G •]. (2.8) 

In practice, we will only be interested in the case that U and X are independent. We will 
sometimes need the elementary relation 

P[Thin^(Thin„(x)) G • ] = P[Thin^„(x) G • ] (v, u G [0, 1]"^, x G N^). (2.9) 

The next lemma gives suficient conditions for one Lloyd-Sudbury model X^ to be a f-thinning 
of another Lloyd-Sudbury model X^. 

Lemma 2 (Thinnings of Lloyd-Sudbury models) Consider {q,a^,b^,c'',d^,e'')-LSM's 
X^ {k = 0, 1,2). Assume that for k = 1,2, X^ is dual to X^ with the duality function from 
Hl-1\) and duality parameter rj^^ . Then 

P[X^ g . ] = P[Thin^(A:^) G • ] implies P[Xl G • ] = P[Thm,,{X^) G • ] (t > 0). 

(2.10) 
provided that v := {I - r?^°)/(l - 7?^°) G [0, 1]. 

Proof This is more or less |Sudnn| Theorem 2.1] (although the formulation of thinning is a bit 
different there), but since the proof is very short, we give it here. We introduce the notation 

ry^:=JJr?^» (r/ G M, x G N^) (2.11) 



Setting 9 := 1 — rj , we see that the duahty between X and X^ can be cast in the form 
E[{l-e''^)^oXt]=E[{l-e''^)^tX^] (t>0). (2.12) 



We need the relation 

E [(1 - 0)Thine, (x)] ^ ^ ^(^ _ qqY] . (2.13) 

If 6* G [0, 1], this is the relation P [Thing (Thing/ (x)) = 0] = P[T\im.ee' {x) = 0]. It is not hard 
to show that 1)2. 13() holds more generally for all ^ € M. Now if the initial laws of X^ and X"^ 
are related as in 1)2. lUj) . then for t > 

^[(1 _ ^10)X0Thm„(X2)j ^ ^j^^ _ ^10)Thm„(X0x2)j ^ ^^^^ _ ^^lO^XgX^j ^ ^^^^ _ ^20^X0x2^ 
= ^[(1 - e20)X0X„2j ^ ^^(^ _ glO)XOThin„(Xg)] ^ ^[(^ _ 0lO)X{>Xij ^ ^[(^ _ ^lO^XOx^j^ 

(2.14) 
Since this is true for each Xq, it follows that Thin„(X|) and X} are equal in distribution. I 

Remarks I did not check the answers to the following questions: 1. (Perhaps not too difficult) 
In Lemma 121 is it sufficient if X^ is only a formal dual, i.e., is it allowed that some of the rates 
oP , . . . ,e^ are negative? 2. (Perhaps somewhat more tricky) Is there a general formula, in the 
spirit of (|2.6p , that tells us when one Lloyd-Sudbury model is a thinning of another? 

3 Contact-voter models 

In this section, we apply the results from the previous section to mixtures of the contact process 
and the voter model. Let q be as in (|2.3j) and let r, s, ttt. > be constants. By definition, the 
contact-voter process with underlying motion kernel q and parameters r,s,m, shortly the 
{q,r, 8,171) -CVP, is the Markov process X = {Xt)t>o in {0,1}^ with formal generator 



GcYpfix) := Y. q{i,j){{r + s)x{i){l - x{j)){f{x + 5,) - f{x)} 

+r(l - x{i))x{j){f{x - 6j) - f{x)}j (3,1) 

+m ^ x{i){f{x - Si) - f{x)}. 

i 

We may interpret Xt{i) £ {0, 1} as the genetic type of an organism living at time t > at the 
site i £ A, where 1 is the fitter type. An organism living at site i invades a site j with rate 
rq{i,j) if it is of type 0, and with rate (r + s)q{i,j) if it is of type 1. In addition, organisms 
of type 1 mutate with rate m to organisms of type 0. We call r the resampling rate, s the 
selection rate, and m the mutation rate. For r = 0, our process is a contact process with 
infection rate s and recovery rate m, while setting s = ni = yields a voter model. 

Contact-voter processes fall into the class of Lloyd-Sudbury models described in the last 
section. Indeed, setting 

a = 0, b = r + s, c = m, d = r + m, e = (3-2) 



in the generator in ()2.2|) yields the generator in H3.1|) . We may therefore apply Proposition^ 
to find duals of contact- voter models. For a given duality parameter rj, Proposition ^ tells us 
that the Lloyd-Sudbury model with the parameters 

a' = 2r77, b' = j^s, c' = m - (1 + 7^)7, 
d =m + j^s, e = -7, 

where 

l = j^s-r (3.4) 

is dual to the (g, r, s,m)-CVP, provided that the rates a' ,b' ,c' ,d' ,e' are all nonnegative (and, 
hence, the LSM with these rates is well-defined). The conditions a',e' > are equivalent to 
T] E (—00,0] U {r/{r + s)}, where r/{r + s) is the value of rj that yields a self-duality. The 
condition 6' > is now trivially fulfilled, while d' > is equivalent to —m/{s — m) < rj ii 
m < s. The condition c' > leads to somewhat messy conditions on ry, but if —1 < ij, this 
is also trivially fulfilled. In the latter case, we can moreover give a nice interpretation to the 
duals with r] < 0. 

Let q be as before and let e £ [0, 1] and p, (3,6 > he constants. Let Y = {Yt)t>o be a 
Markov process in {0, 1} , defined by the formal generator 

GKwfiy) ■■= J2 ^(i^3){py{i){{i - y{j)){f{v + Sj - 5.) - f{y)} 

+yij){efiy - 6, - 6^) + (1 - e)fiy - 5,) - f{y)})) 

+/3y(i) ((1 - y{j)){f{y + 6,) - f{y)} + ey{j){f{y - 5,) - f{y)})) } 

(3.5) 
We may describe the process Y in words as follows. Particles jump from a site z to a site 
j with rate pq{i,j), a particle at the site i gives with rate Pq{i,j) birth to a particle at 
the site j, and particles die with rate 5. If in this process, two particles land on the same 
site, then they annihilate with probability e, and coalesce with the remaining probability. 
We call Y a system of random walks with annihilation, branching, and coalescence, with 
underlying motion kernel q and parameters e,p,P,6, or shortly the {q,£, p, I3,6)-RW. This 
type of models falls into the class of Lloyd-Sudbury models. Indeed, the {q,e, p,P,S)-BW is 
the {q, 2pe, (3, p(l - e) + pe + 5, 6)-LSM. 

Reinterpreting the duals in 1)3. 3(1 . we find the following result. 

Proposition 3 (Duals of contact-voter processes) 

(a) Each (q,r,s,m)-CVP with s > is self-dual with duality parameter r/{r + s). 

(b) A {q,r,s,m)-CVP is dual, with duality parameter —e, to the {q,e,r-\-j^s, j^s,m—j^s)- 
RW, provided that m > jr^s. Conversely, a {q,e, p,P,S)-RW is dual, with duality parameter 
T] = —e, to the {q, p — ef3, (1 + e)/3, 5 + ef3)-CVP, provided that p > efi. 



Note that the formulas in part (b) simphfy a lot if e = 0, i.e., if there is no annihilation. 
By Lemnia|2 it follows that (g, e, p, /?, (5)-RW's are thinnings of each other, and of CVP's. 

Lemma 4 (Thinnings of contact- voter processes) Fix q as in i2.^) . r,m > 0, s > 0, 

and < e < e' < 1 such that m > j^—r^- Then 

(a) The {q, e,r + j^s, j^s, m — ■yj^s)-RW is a (1 + e)~^(l + ^)~^ -thinning of the {q, r, s, m)- 
CVP. 

(b) The {q,s' ,r + i+p'S, i^s,m — jt:^s)-RW is a {1 + e)/{l + e')-thinning of the {q,e,r + 
j^,s,j^s,m-j^s)-RW. 

In particular, it follows that each {q,£,r + jr— 5,^^^5,771, — y^s)-RW is a thinning of the 
(g, 0, r, s, m)-RW, i.e., a system of random walks with branching, coalescence, and deaths 
(with no annihilation). 

4 Local mean field limits 

4.1 Basic definitions 

In this section, we will argue that two types of interacting SDE's, and one particle system, all 
three with applications in population dynamics, can be obtained as 'local mean field limits' of 
contact voter processes and their duals. Moreover, we will investigate how the dualities and 
thinning relations of Proposition |31 and Lemma 0] behave under these limits. Since our main 
interest is in dualities and thinnings, we do not give detailed proofs of our limit relations, as 
this would get too technical, in particular for infinite initial states. Instead, we give the main 
calculations, which for finite A are almost a proof, and roughly indicate what needs to be done 
to make this precise. 

In order to properly define our systems and their duals, we need a few definitions. Let q 
be as in 1)2. 3|) . It is possible to choose strictly positive constants {'Wi)i(z\ in such a way that 

/ ^Wj < cxD and 2.Q{hJ)'Wj ^ Kwi Vi G A (4-1) 

» 3 

for some K < oo. We fix such constants from now on and define a Liggett-Spitzer norm and 
Li-norm 

\\x\\r,:='^Wi\x{i)\, |x|:=^|x(i)| (x e M^). (4.2) 

i i 

Our processes of interest will be Markov processes taking values in the spaces 

£{S) := {xe S^ : ||x|U < oo} and ffin(5) := {x £ S^ : \x\ < oo}, (4.3) 

where S' = N, [0, 1], or [0,oo). 



First, we will consider a stepping stone model with underlying motion kernel q and parame- 
ters r, s, m, shortly the (q, r, s, m)-SSM. This is the Markov process X in <?([0, 1]) or <Sfiii([0, 1]) 
defined by the system of stochastic differential equations (SDE) 

dXtii) = y2q{J,i)iMJ) - Mi)) dt + sXt(.i){l - ^t(O) dt - mXt{i) dt 

(4.4) 
+y/2rXt{i){l-Xt{i))dBtii). 

This model and generalizations have been considered in |SU86j . The parameters r,s,m can 
be interpreted as resampling, selection, and mutation rates. 

Second, we will consider a branching particle system with additional annihilation and 
coalescence, with underlying motion kernel q and parameters a, b, c, d, shortly the (g, a, b, c, d)- 
BPS. This is the Markov process y in <?(N) or £^fin(N) defined by the formal generator 

GsRw/(y) := Y. 9(^' jX^){/(y + ^0 - ^^) - fiv)} 

+a Y, y{i) (y(^) -l)){f{y-25,)- /(y)} + b ^ y(0{/(y + ^i) - fiv)} (4.5) 

+cY,y{i)(y{i) - i){/(y - ^i) - /(y)} + ^J] y(^){/(y - 5.) - /(y)}. 

i i 

This model, with annihilation rate a = 0, has been considered in |AS05j . 

Third, we will consider a super random walk with quadratic killing, with underlying motion 
kernel q and parameters a, /?, 7, where a, 7 > 0, /5 G M, shortly the (g, a, /?, 7)-5'i?VF. This is the 
Markov process Z in £^([0,cxd)) or i?fin([0, cxd)) defined by the system of stochastic differential 
equations (SDE) 

dZtii) = Yq{j,i){Zt{j) - Zt{i))dt + ^aZt{i)dBt{i) + pZtii) dt - -fZt{i)^ dt. (4.6) 
j 

This model has been considered in |HWS05J . 

We will approximate these processes with contact-voter models and their duals, living on 
a lattice of the form A^^ := A x {1, . . . , N}, with underlying motion kernel 

q^^\ii,k),ij,l)) :={ /\ (4.7) 

[^i^ ifi = j,k^l. 

We interpret the collection of sites {{i,k) : k = 1, . . . , N} as a colony of N organisms. Observe 
that 1)4. 7|) says that our underlying motion jumps with rate u^ ' to a site in a different colony, 
and with rate 1 — v^'^' to a different site in the same colony. We will be interested in the case 
that iV — > 00 and v^ ' — > 0, i.e., we have large colonies and most of the interaction takes place 
between organisms living in the same colony. 

When considering dualities, it will be useful to assume that the laws of our approximating 
{0, 1} -valued processes X^ ' are symmetric with respect to permutations of the individuals 
within a colony, i.e., 

P[n,.(xf )) G •] = P[xf ) G •] (i > 0), (4.8) 



for each i G A and for each permutation vr of {1, ... , A^}, where 

(n..(:.o)(i,*):={:|j-;f» 2:L. (") 

is an operator that permutes the individuals in colony i. It is easy to see that if 1)4. 8() holds 
at i = 0, then also for all t > 0. 

4.2 Stepping stone models 

We claim that stepping stone models and branching particle systems can be obtained as local 
mean field limits of contact voter processes and their duals. 

Claim 5 (Stepping stone local mean field limit) Let X^^^ be the {q^^\A^\s,m)-CVP 
and let y(^) be the {q^^\e,A^^ + ji-^s, j^s,m — jj^s)-RW, where q^^^ is as in J^.Tj ^ and 

^(N) ._ i/^rN), r(^) := rN. (4.10) 

for some e £ [0, 1], r > 0, and s,m > 0. Set 

n'=' (4-11) 



Yri^)■.=Eyri^,k), 



fe=l 



let X be the {q, r, s, m)-SSM, and let y be the {q, er, jr^s, (1 — e)r, m — ■^§-)-BPS. Then 
P[X[,^^E-] =^ P[Xq(^-\ implies P[X|^^G-] =^ P[Xt (^ ■] (i > 0), (4.12) 

and 

P[FS,^^G-] =^ P[yoe-] implies P[Y[^^e-] =^ P[yte-] (i>0). (4.13) 

The main ingredients for a proof of Claim [S] can be found in Appendix IB. II The convergence 
in 1)4. 12(1 is quite natural, and shows that the stepping stone model is a good model for gene 
frequencies of an organism living in large colonies, with small migration between these colonies 
(compare the discussion in |SU86I Section 1]). 

We now investigate what happens to the dualities and thinning relations from Proposition|31 
and Lemma in the local mean field limit from Claim [51 For e = 0, the dualities below can 
be found in |AS05I Theorem 1]; for part (b), see also |SU 86. Lemma 2.1]. The case e > in 
part (b) was discovered during work in progress of the present author with Siva Athreya. 

Proposition 6 (Duals of stepping stone models) Let q be as in \2.S\) . r,m > 0, s > 0, 
and e E [0, 1] such that m > jr— •§. Then 



(a) The {q,r,s,m)-SSM is self-dual with duality function 

i;{x, x') ■=e--s^i x{i)x'{i) ^^^ ^/ ^ ^^qq^ ^^y^_ 



(4.14) 



(b) The {q,r, s,m)-SSM is dual to the {q,er, jr^s, (1— e) 



r,m- 



l+S' 



-BPS, with duality function 



il^{x,y):=\{{l-{l + e)x{i)) 



yii) 



(xGf([0,l]), yGf(N)). 



(4.15) 



where we assume that either \x\ < oo or \y\ < oo in case e = 1. 



Proof These dualities can be verified by direct calculation, but the point of this paper is to 
show that they occur as the limits of the relations in Proposition |3J so we follow that road. 

We start with self-duality. Let X and X' be independent (g, r, s,m)-SSM's. We approxi- 
mate X and X' with independent (g'^^r^^^ s,m)-CVP's X^'^' and X (^', respectively, as in 
JUni). By Proposition 01 (a). 



N 



E 



[nn 

i k=l 



-W + , 



.{N) 



rW .Xr>{i,k)X,'''>{i,k) 



'(N), 



E 



[nn 

i fc=l 



.(N) 



'(N), 



rW .Xl,'"{i,k)Xr'{i,k) 



(4.16) 



.(JV) + 



(t > 0). 



Using 1)4. 10(1 . e rewrite the left-hand side of 1)4. 16() as 

,_,^^-NN-^El^Xr{^,k)X^\^,k) 



E 



m 



l + N- 



(4.17) 



Let us assume that the laws of X^ ' and X ^ ' are symmetric with respect to permutations 
of the indivii 
follows that 



of the individuals within colonies. Then, since Xj: and Xq are moreover independent, it 



N 



'(N), 



N-^Y.^ri^,k)x'ri^,k), 



(4.18) 



fc=i 



which is the fraction of the population in colony i that belongs both to Xj: and to Xq , 
converges as A^ — > oo in probability to Xt{i)XQ{i). Therefore, we see that the expression in 
(|4.17|) converges, as A^ ^ oo, to 



[Yle-fMi}^:,ii)], 



E\ Y\ e r 

i 



(4.19) 



The same arguments apply to the right-hand side of (|4.16jl . so we find that X is self-dual with 
the duality function from ()4.14|) . 

To prove also part (b) of the proposition, let A" be a (g, r, s,m)-SSM and let 3^ be an 



independent {q, er, jr^s, (1 — e)r. 



m 



l+e 



)-BPS. We approximate X and y with independent 
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(g W , r W , 5^ m)-CVP's X(^) and (g(^) , e, r(^) + j^s, ^s, m- ^s)-RW's y(^) , as in KTlh 
and KWi . By Proposition El (b), 



^ 



[n 



,J:t,Xri^,k)Yi^\^,k) 



E 



[m-^y 



^N yW 



(N), 



ELlXr\^,k)Y,^"\^,k) 



(4.20) 



(t > 0). 



Assuming that the laws of X^^> and Y^^' are symmetric with respect to permutations of 
the individuals within a colony, we see that Ylk=i-^t (^)^)^o ih^), which is the number 
of organisms in colony i that belongs both to X^ and to Ig' , converges in probability as 
A^ ^ oo to Thmxt{yo)ii)- Treating the right-hand side of (|4.20j) in the same way and taking 
the limit N ^ oo, we find that 



E 



(1-(1 + ^))) 



Thin;t,(3^o) 



E 



(l-(l + e))) 



Thmxo {yt 



(4.21) 



By 1)2. 13() . this means that X and y are dual with the duality function from 1)4. 15|) . I 

There is also an analogon of Lemma 0] for stepping stone models and branching particle 
systems. We adopt the following notation. For y G [0, oo) , we write Pois(y) to denote a 
Poisson measure with intensity y, i.e., a N -valued random variable whose components are 
independent Poisson distributed random variables with mean y{i). If Y is random, then we 
use the symbol Pois(y) for any random variable with the law 



P[Pois(y) £■]:= I P[Y G dy]P[Pois(y) G •]. 



(4.22) 



Lemma 7 (Poissonizations of stepping stone models) Let q be as in \2.t!\) . r,m > 0, 

s > 0, and < e < e' < 1 such that m > ji-rS- Then 

(a) The {q,er, j^s,{l — e)r,'m — j^)-BPS y is a Poissonization of the {q,r,s,m)-SSM X, 
weighted with {1 + e)^^-, i.e., 

P[yo G •] =P[Pois((l + e)-ifA'o) G •] implies P[yt G ■] = P[Pois{il + e)~^^^Xt) G •]. 

(4.23) 

(b) The {q,e'r, j^s, {l—e')r,'m—j^)-BPS is a (l+e) / {1+s') -thinning of the {q,er, j^s, (1— 
e)r,m-j^)-BPS. 

This can be proved by taking the local mean field limits from Claim [S] in the thinning relations 
of Lemma HI Alternatively, these relations can be proved directly using the dualities from 
Proposition ini Since this is very straightforward, we skip the proof. 
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4.3 Super random walks 

We claim that super random walks with quadratic killing can be obtained as local mean field 
limits of contact processes. (Note that in this case, we do not need additional voter model 
dynamics.) 

Claim 8 (Super- RW local mean field limit) Let a, 7 > 0, /? e M, and let X(^) be the 

(gW,0,s(^),m(^))-(7^P, where g(^) is as in Ul^ , and 



sW:=y^^, z^(^):=lA(^\ m(^):=s(^)-/3. (4.24) 

Set (f) := yjapf and 

xf\i^:=-^Y.^'"\^,k). (4.25) 



fc=i 
Let Z be the {q,a, (3,'y)-SRW. Then 

PiX^) e • ] =^ P[^o° e • ] ^rnpHes P[Zf ^ G • ] =^ P[Zf G • ] (t > 0). (4.26) 

N—*co N~*oo 

The main idea of the proof of Claim |H1 can be found in Appendix IB. 21 

The next proposition shows what happens to the well-known self-duality of the contact 
process in the local mean field limit from Claim |S1 

Proposition 9 (Self-duality of SRW) The {q,a, l3,'j)-SRW is self-dual with duality func- 
tion 

^{z,z'):=e-i^i<'>'^'y (4.27) 

Proof We will show that this duality is the local mean field limit of the contact process 
self-duality. Let Z and Z' be independent (g, a, f3, 7)-SRW's. Approximate Z and Z' with 
independent (gW,0, s(^),m(^))-CVP's X(^) and X'(^) as in Claim|Hl By Proposition El (a), 
X^ ' and X ^ ' are dual with duality parameter 0, i.e., with duality function 



^(x, x') = n n 0"^^' '^"'^^' '^ = n \el. ^i.,^wm-o, > (4-28) 



I.e., 



[n 



E \\1 



i fc=l 



^JV y{N)., ,^v'W 






(4.29) 



Assuming that the laws of X^^' and X (^' are symmetric with respect to permutations of the 
individuals within colonies, we observe that 



^x|,''\^,k)X,'^''\i,k) (4.30) 



fc=i 
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converges in probability to Pois(^2^o^t)(0) since for large N, the first process contains approx- 
imately 4 /^-Zo(0-^^ individuals, each having an approximate probability A/^Z^{i)N~'^''^ to 
belong to the dual process as well. It follows that the left-hand side of ()4.29|) converges, as 

A^ — > oo, to 

P[Po[s{lZoZ[) = 0] = E[e-aEiMi)2t{i)y (4.31) 

Treating the right-hand side of ()4.29|) in the same way, we find that Z is self-dual with the 
duality function from H4.14() . I 

We have just seen that Proposition El follows by taking the local mean field limit of Proposi- 
tion |31 (a). One may wonder if taking the same limit in Proposition |31 (b), one may discover 
other duals of super random walks. It turns out that this is not the case. Rather, the dual 
processes from Proposition |31 (b) also converge to super random walks, and one finds duality 
relations that, using scaling, can all be reduced to the self-duality in Proposition EJ Likewise, 
the thinning relations from Lemma|l]converge, under the local mean field limit from Claim|HJ 
to trivial scaling relations. Since this is quite lengthy and does not yield anything new, we 
omit the details. 

A Lloyd-Sudbury dualities 

In this appendix, we outline the proof of a generalization of Proposition ^ Generalizing (|2.2)1 . 
we look at models with generators of the form 

GLsf{x) := 5] {ia(i, j)x«x(j){/(x -6,- 6,) - /(x)} 

'^' +b{i,j)x{i){l - ximfix + 5j) - fix)} 

+c{hJMi)xiJ){fix - 6j) - fix)} (A.l) 

+dii,j)xii)il - xij)){fix - 5i) - fix)} 

+e(i, j>(i)(l - xij)){fix -5, + 6,) - fix)}], 
where a, 6, c, d, e : A x A ^ M are nonnegative functions such that 

sup'^giij) <oo, sup^5((i,j) < 00 ig = a,b,c,d,e). (A.2) 

i ■ j ■ 

Without loss of generality, we may assume (as we do from now on) that aii,j) = a(j, i). 
Generalizing our earlier definition, within this appendix, let us call the Markov process with 
generator Gls the (a, b, c, d, e)-LSM. We will outline a proof of the following generalization of 
Proposition ^ 
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Proposition 10 (Dualities between Lloyd-Sudbury models) The (a, b, c, d, e)-LSM and 
{a' ,b' , c' , d' , e')-LSM are dual with duality parameter rj £ M\{1} provided that 

b'{i,j)=b{j,i) +j{i,j), 

(^{i,j)=c{i,j) -7{j,i) -VlihJ) + {e{i,j) -e{j,i)) +v{KhJ) -^(j,0)> (A.3) 

d'{ij)=d{i,j) +j{i,j) + {e{i,j) -e{j,i)), 

e'{i,j) = e{j,i) -j{i,j), 

where 

l{hj) ■= {a{i,j) +c{j,i) -d{i,j) +r]b{j,i) -e{i,j) +e{j,i))/{l - r]) (A.4) 

Observe that these formulas simphfy to H2.6|) if the functions a, . . . ,e are symmetric. We note 
that an equivalent set of equations for (|A.3|) is 

a'{i,j) +'n{e'{i,j) + e'{j,i)) =a{i,j) +r]{e{i,j) + e{j,i)), 
b'{i,j) + e'{i,j)=b{j,i) +e{j,i), 

i{i,j) = -lij,i), (A.5) 

d'{i,j) + e'{i,j)=d{i,j) + e{i,j), 
e'{i,j) - \i{j,i)=e{j,i) - \l{i,j), 

where ^'{i,j) is defined as in ()A.4|) . with the rates a, . . . , e replaced by a', . • • > e'. 
Proof of Proposition [TUl Bv standard theory |EK86[ Section 4.4], we need to check that 
GLsV',(-,y)(^)=G'LsV',(x, •)(?/) (x,ye{0,l}^), (A.6) 

where G^g is defined as in HA.1|) . with the rates a, ...,e replaced by a',...,e'. Diligent 
calculation reveals that 

= iv'^ - 1)IZ IK^"^"^ " ^)"(''-?') + (^ - V){e{i,j) + eU,i)))x{i)x{j)y{i)y{j) 
+ (a{i,j) -d{i,j) -e{i,j) +r]b{j,i) +c{j,i) + r]e{j,i))x{i)x{j)y{i) 



(A.7) 
l-r])e{t,j)x{i)y{t)y{j) 



+ [d{i,j) + e{i,j)jx{i)y{i) 
-v[KhJ) + e{i,j)jx{i)y{j)j7pr^{x,y). 
and 

G'LS^r, (2;, •)(?/) 

= (7?-^ - 1)^ {^ ((??"' - l)aii,j) + (1 - 7]){e'{i,j) + e'(j,i)))x(i)x(i)y(i)y(j) 

-(1 - r])e{i,j)x{i)x{j)y{i) 

/ / , / , / , \ ^^• 

+ {a'ii,j) -d'{i,j) -e'{i,j) +r]b'{j,i) +c'{j,i) + T]e'{j,i)jx{i)y{i)y{j) 



+ \d'{i,j) +e{i,j)jx{i)y{i) 
-v{b\j,i) + e'{j,i)jx{i)y{j)JTprj{x,y). 
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Comparing (|A.7|) and HA.8|) . we see that HA.6|) is satisfied provided that 

a'ihj) +v{e{'i,j) +e'{j,i)) =a{i,j) +'n{e{i,j) +e{j,i)), 

-(1 -r])e{i,j) =a{i,j) -d{i,j) -e{i,j) +r]b{j,i) +c{j,i) +r]e{j,i), 

a{i,j) -d'{i,j) -e'{i,j) +r]h'{j,i) + c{j,i) +r]e'{j,i) = -{I - r])e{i,j), (A.9) 

d'{i,j) +e{i,j) =d{i,j) + e{i,j), 

b'{j,i) + e{j,i) =b{i,j) + e{i,j). 

which can be solved to yield (|A.3|) . I 

Remark The original motivation for the calculations in this appendix was to generalize the 
self-duality in Proposition IHl (a) to the case that q{i,j) ^ qij,i)- Surprisingly, this does not 
work, unless r = 0. In particular, li m = but r,s > 0, one can check from HA.3|) that a 
{q, r, s, 0)-CVP can only be dual to another {q' , r' , s' , 0)-CVP if q = q' and q{i, j) = q{j, i). 

B Local mean field limits 

B.l Contact voter processes 

Idea of proof of Claim |^ The process X jumps as 

x^x + j^6i with rate ly'^^^N'^A^^ + ^) Yl 90''^)^(i)^(^ " x{i)N) 

+(1 - zy(^))(iV - l)-^(r(^) + s)x{i)N{N - x{i)N) 
x^x-j^6i with rate i^'^^'>N-^A'^^ ^ q{j,i){N - x{j)N)x{i)N (^-1) 

+(1 - zy(^))(iV - 1)-V(^)(A^ - x{i)N)x{i)N 
+mx{i)N. 

It follows that the process started in Xq = x satisfies 

limt-ii?[Xf )« -x«] = ^^(^V(^) Y, qiJ,i)ix{j)-x{z)) 

+ (1 - z^(^))s(l - N~^)-^x{i){l - x(i)) 

—mx{i) (B-2) 

limi^i^[(xj^^(i) - x{i)f] = zy(^)(2r(^) + s)N-'^ ^ q{j, i)x{j){l - x{i)) 

+(1 - z^(^))(2r(^) + s){N - l)-^x(i)(l - x(i)) 
+mx{i)N~ . 
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Under the assumptions H4.1U() . this simphfies to 

hmt-^E[xf\i) - x{i)] = Y, qiJ,^)i<J) - x{i)) 

+sx{i){l - x{i)) - mx{i) + 0(iV"^), (B.3) 

lim t-^E[(xf\i) - x{i)f] = 2rx{i){l - x{i)) + 0{N-^) 

as N ^ oo. This, plus an estimate on a higher moment, can be used to show that the Umiting 
process X started in A'o = x satisfies 

limt-\E[f{Xt)] - fix)) =Y,q{J,i){x{j) - x{i))^^f{x) + s^xi^^^^ 

-m ^ x{i)Q§^f{x) + r ^ x{i){l - x{i))^^f{x) 

i i 

(B.4) 
for all twice continuously differentiable / depending on finitely many coordinates. In the 
right-hand side of ()B.4|) we reckognize the generator of the diffusion process in H4.4|) . More 

formally, the calculations above can be used to show that the processes X are tight, and 



each weak limit point satisfies the martingale problem for the operator in (|B.4|) . Standard 
theory then yields the convergence in (|4.12|1 . 
To prove the convergence in (|4.13|) . set 



(iV)— ^W_^_l_g 

1 

l+e 



P 

P 

6 



h', (B.5) 



■m-j^^s, 



and observe that the assumptions (|4.10j) imply that 






(B.6) 
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ine process Y jumps as 

y^y-6i + 6j with rate u^''^N-'p'^''^q{i,j)y{i){N-y{j)), 
y^y-6i-6j with rate eu^^'> N~^p^^'^q{i,j)y{i)y{j), 
y^y + 6i with rate u^^^N~^P ^ q{j, i)y{j){N - y{i)) 

+(1 - i^WkV- ir'Py{i){N - y{i)), 
y^y-6. With rate (1 - e)i.(^)iV-V^) J] (?(i, j)y«y(i) ^g^^ 

+(1 _ ^)(i _ ^W)(N- i)-V(^)y«(y« - 1) 

+e{l-u^''^)iN-l)-'py{i)iy{i)-l) 
+6y{i), 

y^y-25i with rate e{l - u^^'^){N - l)'^ p^^'^y{i){y{i) - I). 

In the hmit N ^ oo, this simpUfies to 

y ^ y- 5i + 5j with rate q{i,j)yii), 

y ^ y + Si with rate j^sy{i), 

y^y-6i with rate {I - e)ry{{){y{{) - I) + {m - j^s)y{i), 

y^y — 25i with rate ery{i){y{i) — 1), 

plus terms of order N~^ . I 

B.2 Super random walks 

Idea of proof of Claim |^ The proof of (|4.26p is very similar to the proof of H4.12() . Indeed, 
the process X jumps as 

x^x+ -^6i with rate i/^^^TV-^s^^) ^ q{j,i)x{j)(p-^^{N - x{i)(t)'^VN) 

+(1 - z^(^))(iV - l)-h^^'^x{i)^~'^VNiN - x{i)(p-'^^/N) 



(B.8) 



X — > X f^^i with rate m^ 'x{i)4> \N. 



(B.9) 
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It follows that the process started in Xq = x satisfies 

limt-ii?[Xf^(f)-:E(i)]=z.WsW Y^ qiJ,i)xiJ){l-xii)r'N-^/^) 

+ (1 _ A))s(^){xii) - r^N-'/^xiif)il - N-')-' 
—mS 'x{i) 

limt-i^[(Xl^^(i) -x(i))2]=0iV-i/2^WcjW Y^ q{j,i)x{j){l-x{i)^-^N-^/^) 

+0Ar-i/2(i _ i.W)sV)(^(i) _ ^-^N-'/^x{if){l - N-')-' 

(B.IO) 

By H4.24() . for large A'' this simplifies to 



limt-^E[x[''\i) - x(i)] = ^ q{j,i)xij) + (/? - l)x(i) - ^x{i) 



hmt-^E[{xi''\i)-x{i)f]=ax{i), 



- Yl q{J,i)ix{j)-x{i))+px{i)-jx{i)^, (B.ll) 



plus terms of order N ^/^. (Here we have used that ^^ •. j^iq{j,i) = !•) It follows that X 
converges to the process with generator 



{N) 



Gf{z) = Y qU, i)i4J) - 4i))dSr) +l"Yl ^^^J^JTt:^ 



'dz{i) ' 2"Z^-vv5^ 

j^j i 



(B.12) 



i i 

i.e., the process in 1)4. 6p . I 
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